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Abstract 



We present a new method by which wavefunctions with simple structure 



^ \ are renormahzed so as to contain more comphcated structure. This method, 

H ; 

5^ , caUed ff^-cooUng method, is apphed to the study of the quadrupole collective 

motion of ^^Fe, ^^Cr, ^^Fe and ^^Cr. The shell-model wavefunctions of lowest- 
lying states of these nuclei are well treated by this method. By using the 
wavefunctions obtained via the i?"-cooling method, IBM-2 parameters are 
derived from a realistic shell-model hamiltonian and transition operators. The 
Majorana interaction becomes sizably repulsive, primarily as an effect of the 
renormalization. The bosonic E2 effective-charges are enhanced due to the 
renormalization, while a quenching occurs in the Ml and MS parameters for 
proton bosons. It is shown that the x parameters take similar values in the 
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hamiltonian and in the E2 operator. 
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I. INTRODUCTION 



Middle pf-sheW nuclei provide us with a precious testing ground to understand various 
aspects of the quadrupole collective motion from microscopic standpoints. Computational 
difficulties in a realistic shell model rise rapidly in general, as the mass number increases. The 
growing computer power, however, enables us to carry out realistic shell-model calculations 
in the middle p_^shell region. On the other side, the middle p_^shell nuclei seem to gain 
significant quadrupole collectivity, which is a global and dominating feature of heavier nuclei. 

Recently Sebe and the authors have reported one of the most successful shell-model 
results for = 28 — 30 nuclei The Kuo-Brown interaction [Q, which had been derived 

from a realistic A^A^ potential through G-matrix, has been employed in these calculations, 
together with a large configuration space including excitations from the O/7/2 orbit. To be 
more precise, considering the following configuration, 

where ni = {Z - 20) + 8 and r^a = - 28 for the 20 < Z < 28 < A^ < 40 nuclei, we have 
adopted a space consisting of all the k = 0,1 and 2 configurations. It has been confirmed 
that, for even-even nuclei, the energy levels are reproduced remarkably well for < 4MeV. 

The presence of mixed-symmetry states with respect to the proton and neutron collective 
degrees of freedom has been predicted by the proton-neutron interacting boson model (IBM- 
2) It has been pointed out P] that a mixed- symmetry 2+ state may lie lower than the 
other mixed-symmetry states in spherical nuclei, although the origin of such a low-lying 
mixed-symmetry 2+ level has remained open. Experimental studies have suggested that the 
mixed-symmetry state exists around = 3MeV in the Cr-Fe region 0. A realistic 
shell-model analysis has been applied to pin down the mixed-symmetry 2+ state of ^^Fe [0 , 
clarifying which states share substantial fractions of the mixed-symmetry component. It is 
of special interest to study this type of collective modes more extensively, on the basis of a 
realistic shell-model calculation. 



3 



The realistic shell-model hamiltonian couples collective degrees of freedom to non- 
collective ones, in general. When large-scale shell-model results are interpreted in terms 
of lBM-2, it is important to incorporate, through a certain renormalization procedure, ef- 
fects of relevant non-collective degrees of freedom into the calculations made in the collective 
subspace. This is an example of the general problem as to how a complicated system can 
be described with a limited number of degrees of freedom by taking into account a variety 
of correlations in an effective manner. Rayleigh-Schrodinger's perturbation theory con- 
stitutes a possible way, by which the model wavefunction is modified. The second-order 
perturbation has been applied to renormalize the lBM-2 parameters [§,0. Another way is 
Bloch-Horowitz's renormalization of operators, in which operators, rather than wavefunc- 
tions, are modified perturbatively so as to carry the relevant correlation effects. Some useful 
general theories have been developed by extending Bloch-Horowitz's method: Feshbach's 
projection method [jlO[ and the folded- diagram theory for example. These methods are. 



however, more or less based on the perturbation theory. In the cases to be considered in 
this paper, perturbative ways are inappropriate, as is argued just below. 

Our present goal is an investigation of the quadrupole collective states which are to be 
described within IBM-2, in connection with the realistic shell model. In the first approxi- 
mation, the s- and d-bosons in IBM-2 correspond to the collective 0"'' (S) and 2+ (D) pairs 
of valence hke-nucleons In Cr-Fe nuclei, as is assumed in Ref. 0, the S- and D-pairs 
normally comprise only the k = configuration of Eq.(|l]), because of the Z = N = 28 
magic number. The realistic shell-model wavefunctions, however, contain other configura- 
tions. According to the reahstic shell-model results, the leakage out of the A; = space is 
so significant that even the Oi" and 2i wavefunctions are not well enough covered with this 
usual SD space (< 60%) 0. In order that the Of and 2i states can be described within the 
IBM-2, the correlations beyond the SD-pairs must be taken into account. This relatively 
large A; > fraction prevents perturbative ways from being applicable. A method beyond the 
perturbation theory is required. It is commented that this situation takes place because the 
^^Ni core is not very stiff. Perturbative approaches may be legitimate in other mass-regions 



in connecting IBM-2 to realistic shell-model. 

We recall here that, as far as several lowest-lying levels are concerned, they are success- 
fully reproduced by the Horie-Ogawa hamiltonian ||12| with only the k = configuration 
in the Cr-Fe region, apart from the precise description of the mixed-symmetry states [|I|. 
Moreover, this k = shell-model result is connected with IBM-2 fairly well, at least for O^*" 
and 2i [|l3l. This fact suggests that, even in more realistic cases with A; > configurations, 
the lowest-lying states may be described within IBM-2 through a proper renormalization. 
In Ref. [H] modified S'D-pairs have been introduced, and the fragmentation of the mixed- 
symmetry 2+ components was clarified for ^^Fe. 

We introduce a new and yet simple method in this article. It is applicable even to some 
cases where the perturbation does not work well. The method is applied to quadrupole 
collective states of Cr-Fe nuclei, and the wavefunctions of those states are renormalized. In 
addition to ^^Fe and ^"^Cr, on which the shell-model results have already been reported in Ref. 
0, ^^Fe and ^^Cr are studied. Furthermore, by extending the OAI mapping 0, the IBM-2 
hamiltonian is derived from the realistic shell-model hamiltonian. This is the first work 
of this sort, while there have been many works evaluating IBM-2 parameters from more 
schematic interactions, for instance the surface-delta interaction. The IBM-2 transition 
operators are obtained as well. Renormalization effects on various IBM-2 parameters are 
discussed. Focusing on the IBM-2 results more concisely, we shall investigate properties of 



the mixed-symmetry states in Cr-Fe region in the following paper [14 



II. i?" COOLING METHOD (H"CM) 

The present renormalization method is introduced in a general form, in this section. 
Some details of the procedure will be illustrated in Section |IT|. Although this method may 
be applicable to other many-body problems, we shall apply it, in this paper, to elicit a 
collective space out of the shell-model space. This collective space should correspond to 
that of IBM-2. 
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In ^^Fe, for example, we have a pair of proton holes and a pair of valence neutrons within 
the description assuming the ^^Ni core. The S- and D-pairs of protons are defined as the 
0^ and 2+ states of the (0/7/2)"^ configuration, while those of neutrons are collective 0^ and 
2+ states of the (0/5/21^3/21^1/2)^ configuration. Once the structure of the neutron 5"- and 
i^-pairs is given, the SD space is constructed by these proton and neutron S- and D-pairs. 

We first introduce a subspace of the original Hilbert space. This subspace is denoted 
by W^^\ In the application discussed in Section the original space corresponds to the 



shell-model space, and W^^^ to the SD space. The bases belonging to W^^^ are hereafter 
called primary bases, while those outside W^^^ non-primary bases. It is required that the 
primary bases include basic dynamics already. If there is a conserved quantum number J, 
the space W^'^^ can be decomposed as 

l^(o)=01^f. (2) 
J 

In the practical case, J represents nuclear spin. In the following procedure, wj°^'s with 
different J's never mix with one another, reflecting the conservation law. 

We consider a primary basis \1'^'^^ G Wj'\ where A is the index of basis-state vectors for 
a given J. The quantum number J is not explicitly shown in \E'^°^ for brevity. The state 



"^^^^ evolves with time t as e 



^1°^), where H means the original hamiltonian. When 



A 

the inverse temperature (3 = it with imaginary t is employed, the time evolution e~*^* is 
converted to cooling e~^^. In order to simplify the following discussion, we assume without 
loss of generality that all the eigenenergies are non-negative. This situation is attained, if 



necessary, by shifting the origin of energy. The expectation value (\E'^°'' 



^1°^), which 



is a function of /3, is a superposition of exponentially-decreasing components corresponding 
to eigenvalues of H. The number of these components is much larger, in general, than 
the dimension of W^-^\ because a sizable fraction of e~^^ ^1"'') escapes out of VT^*^) with 
increasing 

p. ^-(3H i^jo) ^ ^jo) therefore e"^-^ ^ W^^l We consider, in this 

paper, the situations in which the primary bases form a major part of some low-lying states 



and then (\1';^ 



(0) 



"^^x^"^ is dominated by one or a few slowly-decreasing components, 
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while fast-decaying components are superposed with far smaller amplitudes. If we choose 
appropriate states from non-primary bases, only a small number of them will have a sizable 
mixing with "^^x^ . These non-primary but relevant bases are here expressed as (p^^^ . The 
superscript (u) is used as an index of degree of the coupling to "^^x^ whose meaning will 
be specified later in this section. By taking into account the influence of the (p-hases, the 
wavefunction of "^^^^ will be renormalized as 

^xoc¥;:^+j:cuA'\ (3) 

V 

where c^^x represents mixing amplitude of the 0-basis. The basis "^x is constructed so as 
to contain higher-energy components with significantly small amplitudes. By doing this, 
rapidly-decreasing components in (\['a| e~^^ |\E'a) can be made negligibly small. Then the 
exact low-lying eigenstates will be reproduced to a good approximation by appropriate linear 
combinations of ^E'a's- The following discussion will exhibit how to choose (^^^ efficiently and 
how to evaluate c^^x- With an adequate set of 0a^'''s and Cj^^a's, the truncation up to relatively 
small V in Eq.(|^) is expected to yield good renormalized bases, as will be shown with concrete 
examples. 

In the present method, it is required that the Hilbert space Wj = {\E'a; A = 1,2, • ■ ■}, 
which consists of the renormalized bases, fulfills the approximate relation 

e-^"" Wj ^ Wj, (4) 

up to a reasonably large p. By defining the total space of the renormalized bases by 

W^^Wj, (5) 
J 

Eq.(||) can be expressed as 

e-^^ W ^W. (6) 

We look for a W which satisfies Eq.(^ and, at the same time, remains a rather small subspace 
of the full Hilbert space. Eq.(D or (|^) indicates approximate closure of the renormalized 
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space. The closure is exactly satisfied if W consists of eigenstates of H. Such a construction, 
however, only means a calculation in the full space. We are here seeking to obtain Wj with 
a limited number of the (jy-hases in Eq.(|^), discarding degrees of freedom coupled to the 
primary bases only weakly. For this purpose, we shall trace how the original basis 
evolves by e~^^ . 

We consider a small (3 by rewriting (3 as A/3 for the time being, although it is not essential 
as discussed later. By expanding e~^^^ into the power series of A(3, the cooling of "^^^^ gives 



— i lA V 



, (7) 



where O {{APY) represents an operator with the order of We define the (f)\ bases 

(see Eq.(^)) from the right-hand side. 



M/f). (8) 



Here stands for an appropriate orthonormalization, whose concrete definition will be 
given in Section |T|. The basis <p^^'^ directly couples to "^f^ via H, exhausting the coupling 



leading out of Wj^\ The next basis (j)^^' affects the primary state via its coupling to 



In this manner, important bases are extracted one after another. Note that all the 0-bases 
carry the same quantum number J as "^^x^ ■ 

Since the shell model is defined as a finite dimensional many-body problem, the cooling 
operator e~^^^ does not need infinite series expansion. Moreover, since we have postulated 
that the primary bases include basic dynamics, the number of relevant degrees of freedom 
which couple to the primary bases can be relatively small. Making good use of these features, 
g-A/j/f handled by the power-series expansion as in Eq.(|^), and the 0-bases are generated 
by Eq.(l). 

It has been known that the Lanczos diagonalization algorithm is efficient to obtain 
eigenenergies and eigenfunctions of low-lying states. We here try to utilize the advantage of 
the Lanczos method. The Lanczos method can be derived via the power-series expansion 
of e~^^^ acting on an arbitrary basis. Besides a difference in the P*^ operator (see Sec- 



tion PI) , the 0-bases in Eq.(|D have basically the same form as the Lanczos bases generated 



from . However, it is a key point of the present method that the primary bases contain 
basic dynamics of the system. In other words, the primary bases form the main part of 
the wavefunctions of the low-lying states under interest. Thereby we can regard the current 
procedure as a renormalization. The relation of the present method to the Lanczos method 
will be discussed further in Section |VI[ 

Let us begin with a simple case in which there is just a sing le state in W^j^l We 
do not need the label A in this case. The operator in Eq.(^ expresses the Gram- 
Schmidt orthogonalization to '^^'^^ and (p^'^'^^s with u' < v. We thus generate a subspace 
r(") = 1^(0), 0(1), 0(2), . . . , 0(")|, corresponding to the order of {1^(3^. 

A renormalized wavefunction (see Eq.(|])) is introduced within r^"-). 



^(") oc ^(0) + Cy(l) 



(9) 



with the amplitude c,^'s to be determined. The submatrix of H for the subspace F^") is 
constructed, and the eigenvector associated with the lowest eigenvalue is adopted as the 
renormalized basis \&(") in Eq.(^. It is noticed that the mixing amplitude Cy thus obtained 
depends on n, though this dependence is not explicitly shown here. The basis \E'(") yields 



-A/3_fi" 



vl>(") 



*("))+ fo((A/3) 



(10) 



where £"(") = (\E'(") H \&(")). By this procedure, rapidly-decreasing components become 



substantially smaller in /vl>(") 



-A/3// 



-A/3// 



^(o)V The ^(°) state i 



IS 



^(")) than in (^(o) 

thus cooled down. Increasing n step by step, we can monitor what components are adopted 
in higher-order steps. 

When we have /(> 1) bases in wf'\ the renormalized basis ^^""^ is obtained for each A by 
diagonalizing a submatrix of H in r^"\ which is generated from ^1^^°^ The \E'^"'' (A = 1, ■ ■ ■ , /) 
bases span a space PVj"''. However, H may produce cross-over couplings among bases in F^"'^ 
with different A, which gives rise to a non-orthogonality between bases with different A. In 
order to avoid this, F^""* is created so as to be orthogonal to F^7^ if A 7^ A', by carrying out 



an orthonormalization. This orthogonalization will be illustrated concretely in Section |T1 



Although this modification can generally break the relation (^0|) for individual basis, the 
orthonormalization can be made (see Section [TTTD so that a similar condition should be 



satisfied for the space W^'^^ 
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Af3H ^(n) ^ ^(n) ^ ((A/?)"+l)] . (11) 

This indicates that VF"]"^ fulfills Eq.(|) up to O ((A/5)"). The entire space of the renormalized 



bases at the order n is then defined as 



H^W=01^j"). (12) 

The cooling Eq.(|lT]) is carried out step by step, through the power-series expansion 
(0). We shall call the present procedure if "-cooling method (H"CM). The H"CM gives a 
wavefunction renormalization, incorporating dynamical correlations contained in H. As far 
as the H^CM is concerned, the state e~^^ ^1°'') with small (3 is decomposed in terms of 
the n = 1 basis and the rest. The latter has higher energy than the former, giving 
rise to the faster-decreasing component. Analogously, for a general n, the H"CM process 
produces n faster-decreasing components in addition to the slowest-decreasing component 
{i.e., "^^x^). The closure of the renormalized subspace W^"^^ is fulfilled up to 0((A/9)"), as 
is shown in Eq.(|TT]). The larger n assures the better approximation from the viewpoint of 
the condition (|) or (|). If converges with n, no new basis is created by H acting on 
this subspace. The convergence then becomes independent of (3, which means that Eq.(^ 
holds for a general value of P, not only for Af3. 

Although a cooling can be made only by operating e~^^^ on "^^^^ H"CM the 

diagonalization is performed within T^^'' each step. This accelerates the cooling to an 
appreciable extent, since the diagonalization is equivalent to the full cooling within the 
relevant subspace. Moreover, as far as the dimension of H is finite, bases are exhausted 
at finite n. Therefore, the P —>■ oo limit, which is required for the full cooling in infinite- 
dimensional cases, is not necessary. Because of these properties, all the major components 
for low-lying states are generated with relatively small /3, and the H^CM is expected to be 
efficient even with rather small n. We shall see it in practice in Sections 
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Here we should add the comment that some nuclear collective states, for which the H"CM 
will be used, are not necessarily the lowest-lying state with a specific spin-parity. In such 
cases the term 'cooling' may not be appropriate, and some caution will be necessary in 
applying the present method. A prescription will be shown in Section |VT|, while the actual 
case of the Cr-Fe nuclei will be presented in Ref. M. 



III. ILLUSTRATION OF H"CM 

The H"CM is illustrated in some detail with an example: let us consider the set spanned 
by the S'D-states with = 2+ in ^^Fe. 

As has been mentioned in the preceding section, the proton S- and D-pairs have the 
(0/7/2)"^ configuration in ^^Fe, while the neutron pairs have the (0/5/21^3/21^1/2)^ config- 
uration. The set of 2+ states within this SD space of ^^Fe comprises the following bases. 



2+{D')-F = l) = [\D^)^\D^)f\ (13b) 

2+{SD); F = 0) = i\D^) ^\S,)- \S^) ® |D.)) . (13c) 

These bases straightforwardly correspond to the IBM-2 bases through the OAI mapping 
p. In the following discussions, the proton-neutron property of the wavefunctions is taken 
into consideration, so that the states should correspond to the IBM-2 states with good 
-F-spin. The variable F on the left-hand side of Eq.(|13]) indicates the F-spin value of the 
corresponding IBM basis after the OAI-mapping 0. The maximum of F is obtained by 
-^max = |A^^) where A^^ = A^^ + is the total number of ^SD-pairs, which is the same 
as the total boson number in the IBM-2, for each nucleus. The states with F = F^ax 
are called totally-symmetric states in the IBM-2, while those with F = Fmax — 1 mixed- 
symmetry states. We refer here the bases in the SD fermion space in an analogous manner. 
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Note that F^ax = 1 in ^^Fe. Though the above \2+{SD); F = 0) state is a totally anti- 
symmetric state, it is called mixed-symmetry state in this article, because it belongs to the 
class of the F = Fmax — 1 states. 

For each nucleus we shall consider the complete set of orthonormal bases belonging to 
Wj^^ {i.e., the SD space with a specific spin); '^2'\ ■ ■ •, '^f'^ ■ In the above case of 

Eq.(|l3D, the ^(°)'s turn out to be 

2+(5D);F = 0). (14) 

The ordering of the bases may affect the process of the H"CM, as will become transparent 
below. In this example, the S'D-bases are ordered so that the class of states with full F-spin 
symmetry (F = Fmax) should come first, those with next highest F-spin (F = Fmax — 1) 
come second, and so forth. Within each sector of a given F-spin, the bases are placed from 
the lower seniority to the higher, similarly to the OAI mapping. 

In the H"CM, the 0-bases of Eq.(D as well as the ^(°)'s are generated in the following 
order. 





*f . ■ 






4". ■ 


■. tf' 






■. W^' 


4"', 


4"', ■ 


4'" 



(15) 



Recall that (p\ is generated from H'^'^\ , apart from the orthonormalization by P . We 
define the P*^ operator in Eq. (||) as follows. The first / bases are 

M°\vl/f, vl/W, which 

are already orthonormal, and P^ acts as the unity for them. The (/ + l)-th basis is cp^^'' 
generated from H'^f'^ with the Gram-Schmidt orthogonalization to ^2'''' 5 ' y "^f^ ■ 

Namely, 

rf^-^f<),.,^(o)J■^^f^ (16) 



^1 



(0) 



2+(5D);F 



2+(F'2);F = 1 
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where P^y represents orthogonahzation to the states specified in the curly bracket, together 
with the normalization. The (/ + 2)-th basis is created similarly, except that it should 
be orthogonal also to (p^i\ 

One can repeat the procedure until all the orthonormal bases in Eq.(^) are obtained. Each 
0-basis can be represented explicitly as 

Note that, as n becomes larger, some bases may vanish due to the orthogonalization. The 
bases in Eq.(|T5|) are produced in this manner, by carrying out the Gram-Schmidt orthogo- 
nalization on them successively. 

We then consider a subset T^"^ = 0^^^ 0^ , ■ ■ ■ , ^JT^} for each A(= 1,2, ■■■,/). 

It should be noticed that F^"^ is spanned by the bases constituting the A-th column of 
Eq.(P^. In order to obtain a renormalized basis {i.e., to evaluate c^^x of Eq.(^), we 
construct a submatrix of H within this subspace r^"\ After diagonalizing this submatrix, 
the lowest eigenstate is taken as the A-th basis in Wj^\ In this manner, is obtained 
from T^^^ = {¥^\(j)^^^}, from T^f^ = [¥^\(f)^^\(f)^^^], and so forth. The H"CM space 
W^'^^ is spanned by the bases thus obtained. 

The H'^CM will be useful for extracting some simple structural features from compli- 
cated shell-model wavefunctions. Since the renormalized wavefunctions of the SD states are 
explicitly constructed, it is possible to compare them directly to the shell-model wavefunc- 
tions. It is also straightforward to evaluate matrix elements of a given operator in the space 
W'-''\ Although the lar ger n implies the better closure of the subspace from the viewpoint 
of Eq.(^), we consider n < 2 cases in the following application. 
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IV. APPLICATION OF H'^CM TO SD SPACE IN Cr Fe NUCLEI 



The H"CM is applied and tested numerically in Cr-Fe nuclei, starting from the S'D-pair 
states. 

As has been shown in Ref. 0, the shell-model calculation with the Kuo-Brown realistic 
hamiltonian in the k < 2 space successfully reproduces the observed states up to i?^. ~ 4MeV 
in ^"^Cr and ^^Fe. While the proton S- and D-pairs are uniquely determined by the (0/7/2)^^ 
configuration, the structure of the neutron pairs has to be fixed. In ^^Fe, the structure of 
\S^) is determined so as to maximize the overlap between the 15"^) (S> \ Siy) state and the shell- 
model Oi state. The structure of \Di,) is determined so that the overlap between \St,) {D,^) 
and the shell-model 2^ state should be maximum. The neutron pair structure for ^^Cr is 
fixed so as for \S^) ® \S^) {\S^) (g) \Du)) to have maximum overlap with the shell- model Of 
{2i) state. The structure of neutron pairs is slightly different between ^^Fe and ^'^Cr. Note 
that the seniority projection [Q is carried out when J) is produced in Cr. 

The H"CM is applied, starting with the collective space composed of these SD-pairs. 
The H"CM is driven by the shell-model hamiltonian in the k < 2 space. The subspace Wj'^ 
is spanned by the products of the above S'D-pairs with angular momentum J. Any basis in 
jy{o) carries the lowest isospin, and the isospin is conserved during the H"CM. The primary 
bases \l/^°'''s in Wf'^ are put in the order of F-spin and seniority, as has been mentioned in 
the preceding section. 

For n > 1, \l/^"'''s imply renormalized SD states, while they are referred to by the 
corresponding ^SD-pair structure of "^f^- Once \E'^"^'s are obtained, eigenstates within 
Wj"^ are calculated by diagonalizing the submatrix of the hamiltonian whose elements are 
In Fig.|l|(a), the energies of lowest O"*" and 2+ states within the subspace 
are shown, in comparison with those of the shell model, for ^^Fe. The origin of energy is 
set to the shell-model ground-state energy. The unrenormalized S'D-states gives far higher 
energies (the most leftward sector) than the shell-model eigenenergies (the most rightward 
sector). In this regard, the bare S'D-states {i.e., the primary bases) are insufficient. The 
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renormalization via the H"CM reduces this discrepancy of energies quite efficiently. The 
H^CM appears to recover the shell-model Of and 2i states satisfactorily well. This ob- 
servation is confirmed by direct comparison of the wavef unctions. Table | shows overlaps 
between the S'D-states and the shell-model eigenstates for the lowest 0"*" and 2"*" states. The 
unrenormalized states are certainly different from the shell-model eigenstates. A large part 
of the discrepancy comes from the k > configurations. On the contrary, the wavefunctions 
after the H^CM are quite close to the corresponding shell- model ones, having more than 
90% overlaps. 

Figure [l|(b) depicts energy levels for which energies are measured from the ground state 
defined in each space. It is remarked that the levels without the renormalization resemble the 
ones after the H^CM, whereas the H^CM spectrum is certainly different. The 2f excitation- 
energy in the original SD space is in good agreement with the shell- model result, and 
therefore with experiments. In the H^CM result, the Of state is greatly lowered owing to 
the coupling to non-SD {i.e., non-primary) degrees of freedom. Though the 2f state is also 
lowered, the non-SD effect is smaller than in 0^. Some additional non-SD effect is absorbed 
by the H^CM, which recovers the 2f excitation energy. In an analogous manner, as far as the 
excitation spectra are concerned, the result after the H^CM is close to the unrenormalized 
one for most lower-lying states. There is a certain difference in the l"*" state. We shall 
return to this point later. The 4f state appears to be too high, even in the H^CM result. 
This state seems to be largely influenced by the G-pair degrees of freedom. Although some 
parts of them are included in the renormalized wavefunctions, they are not yet sufficient for 
compensating the whole influence of the G-pairs. 

In Cr nuclei, there is no \Dl; 0+) basis {i.e., 0"^ state with (0/7/2)^ and seniority 4). It is 
not always possible, thereby, to create \E'^°^ basis having a good F-spin value. For instance, 
the \0^{SD'^)) bases with good F-spins are 

0^{SD'y,F = I) = ^ (1^'; 0+) ® \S,) + V2 [\SM ® \Du)f^) , (19a) 
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0+{SD'y, F = l) = ^{^/2\Dl■, 0+) ® - [\SM ® 1^.)] 



(0) 



(19b) 



Since |D^;0"*") does not exist in the present case, only a single \0^{SD'^)) basis is possible, 
and we introduce 



0+{SD^y,F 



(0) 



(20) 



We replace the higher F-spin basis ( |19a|) by (pOl). The lower F-spin basis is also subject to 
similar changes, and should be modified with proper orthogonalization to the basis assigned 
with higher F value. This orthogonalization, however, annihilates the basis corresponding 
to Eq.( |19b|) , because of the lack of the |D^;0"^) component. The following primary bases 
are thus obtained, 



0+(5D2);F=^^' 



0+(D-^):F= - 



(21) 



The \2^{D^)) bases are handled in an analogous manner. The energy levels in the SD space 
of ^"^Cr thus constructed are shown in Fig.^. 

For the = 32 nuclei, the shell-model energy levels of ^^Fe and ^^Cr are presented 
in Figs.^ and ^, in comparison with the observed ones. This shell-model calculation is 
performed by using the computer code VECSSE [T^. The k < 2 space leads to the M- 
scheme dimension of 631, 670 for ^^Fe and 621,478 for ^^Cr. Because available experimental 
data for these nuclei are not as abundant as for the = 30 isotones, a stringent assessment 
of the calculated results appears to be difficult. However, we see that this shell- model 
calculation reproduces the observed levels reasonably well also for these = 32 nuclei. 

The structure of neutron S- and D-pairs in ^^Fe (^^Cr) is assumed to be the same as 
in ^^Fe (^^Cr), for the sake of simplicity. The \Dl; J) bases are produced so as to have the 
generalized seniority |0 of four, by removing lower- seniority components. In ^^Cr, the lack 
of the \Dl;0'^) component causes a modification of the SD bases, as in ^^Cr. The energy 
levels in ^^Fe and ^^Cr within the SD space are shown in Figs.^ and ^ in comparison with 
the shell-model ones. 
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Fi gures ^ 1^ and ^ indicate, respectively, that the H^CM works well for the 0^ and 2^ 
states of ^^Cr, ^^Fe and ^^Cr, to the same extent as in ^^Fe. As shown in Table |, the 
wavef unctions, as well as the energy levels, are in good agreement with the realistic shell- 
model ones. In ^^Cr and ^^Cr, the H^CM wavefunctions have more than 95% overlap with 
the shell-model eigenstates for 0^, and more than 90% for 2^, as in ^^Fe. These numbers 
are somewhat smaller in ^^Fe, originating in the smaller overlaps of the unrenormalized in- 
states and the shell-model eigenstates. Though we have presumed the same neutron-pair 
structure as in ^^Fe, a different choice may improve the overlaps. The overlaps of the H^CM 
wavefunctions with the shell-model ones are still large, exceeding 90% for O^" and 80% for 
2f. Thus the quadrupole collective motion in these nuclei can plausibly be described via 
the H^CM, based on the realistic shell model. 

Although the calculation without the renormalization yields quite high ground-state 
energy, it produces the energy spectrum similar to that in the H^CM result for most lower- 
lying states. This consequence is not a trivial one, and seems to give a rather deep insight 
upon the correspondence between the iSD-pair states and the IBM-2 states. 

It should be mentioned here that the renormalization proposed by van Egmond and 
AUaart [^] has similar aspects to the present one. They have included various correla- 
tions induced through a nearly realistic shell-model hamiltonian. To obtain a renormalized 
wavefunction of the state with one D-pair, the hamiltonian is diagonalized on the one- and 
two-broken-pair bases. Most of these bases are involved in the H^CM, if their work is com- 
pared with the H"CM. However, a few of them emerge only in H^CM. In this sense, a part 
of the H^CM effect has been taken into account in Ref. [0]. The II"CM exploited in the 



present work, on the other hand, provides us with a systematic way to pick up important 
bases. Wavefunctions of higher SD-states are also handled by the II"CM. A significant 
difference from Ref. [0, as well as from other studies, is that the ^^Ni-core excitation plays 
a certain role in the present case, which may be characteristic to the present mass region. 
We should note that, because of this point, the 0~^{S^^); F = F^ax) basis is modified to an 
appreciable extent. 
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It has been demonstrated how efficient the H^CM is. Despite the relatively heavy leakage 
out of the SD space {W^'^^), the H^CM yields reasonable energies and wavefunctions for the 
lowest-lying states. The leakage mainly arises from the k > configurations, namely from 
the ^^Ni-core breaking. It is commented that, if we had an appropriate effective interaction 
in the k = space, the H^CM, or even the bare SD wavefunctions, might have worked more 
efficiently. 



V. IBM 2 PARAMETERS 

We next calculate lBM-2 parameters by extending the OAI mapping, based on the 
H"CM. There have been numerous investigations devoted to the derivation of the IBM- 
2 hamiltonian from schematic interactions like a pairing-plus-quadrupole interaction or a 
surface-delta interaction P,P|,|2D[] . A semi- realistic interaction with the single-range Yukawa 



form has been applied to investigate several of the IBM-2 parameters in Ref. ||19l. On the 
other hand, this is the first work to derive IBM-2 hamiltonian from a realistic shell-model 
interaction; the Kuo-Brown interaction, in the present case. 

We employ the following form of the IBM-2 hamiltonian, including all the possible one- 
and two-body terms, 

P='^,'^ P=T^,i^ J=l,2,3 

+ hj[didlf ^ ■\dJuf \ (22) 



J=0,2,4 



where 



Qp = [dlsp + sldp] + Xp [dldp] , (23) 
y? = IvoAsWp ■ Idpd.f'' + h.c.} + ^V2Adlsl ■ [d^f^ + h.c.} 

+ 1 E cjAdldlf^ ■ Kd,f\ (24) 

^ J=0,2,4 

The so-called Majorana terms, which control the energy of the mixed-symmetry states rel- 
ative to the symmetric states, are defined by 



Mj = [44]^^^ ■ [44]^^^ (for J = 1, 3), (25a) 

M2 = \[dlsl - 44] ■ [4^^. - sM''^ (25b) 

In the OAI mapping, a boson image of a certain nucleon operator is obtained from matrix 
elements concerning low-seniority states. The parameters are fixed so that a boson matrix 
element should be equal to the corresponding matrix element in the collective fermion space. 
We here consider F-spin also, in addition to the seniority. The boson matrix elements are 
equated to the fermion ones, similarly to the sequence of the bases in the H^CM discussed 



in Section |T|. This procedure is briefly illustrated below. 

We first consider the s-boson condensate |s^^). Since = [s^^\H^\s^^), the parame- 
ter Eq in Eq.(p2D is fixed by the equation 



Eo = (0+(5^'); F = F^,, H 0+(5^"); F = F^,^ , (26) 



where H stands for the shell-model hamiltonian. When we regard the |0'^(S'^'^)) basis as the 
unrenormalized one, we obtain the unrenormalized value of Eq from Eq. (pB]) . By putting 
the renormalized wavefunction for |0+(5'^^)), the renormalized value of Eq is evaluated. 
In this procedure, the wavefunction renormalization for the collective fermion states gives 
rise to a renormalization of the IBM-2 parameter. The n parameter is fixed from the 
{0'^{S^'')\H\0+{S^''-^D^)) matrix element. Then e^, and are determined from the 
following coupled equation, 

1 2 A^^ A^^ 



2+(^^ -^D); F = F^,. H 2+(5^°-iD); F = F^,.) , (27a) 



2+(^^"-iD); F = F^,. H 2+(5^°-iD); F = F^,. - l) . (27b) 

The other parameters are evaluated in an analogous manner. It should be noticed that in 
the unrenormalized case this procedure yields the same results as the OAI mapping. 
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The resultant parameters, the unrenormahzed ones and the renormahzed ones via the 
H^CM, are displayed in Table y. For neither ^^Fe nor ^^Cr, the interaction among neutron 
bosons {Vj^) has any effect, since A*"^ = 1. Hence the parameters in V,^ have not been 
determined for these nuclei. Likewise, in the Fe nuclei is not given. For the Cr nuclei, 
the lack of the \Dl; 0+) component is realized by setting cq^tt = oo, while Vq^t, is indeterminate. 

It is often assumed, in phenomenological studies, that the IBM-2 hamiltonian is com- 
prised only of the second, fourth and fifth terms of Eq.(^). Note that, as far as excitation 
energies are concerned, the constant term Eq plays no role. The 2f excitation energy is 
governed mainly by the Nd and Qtt ■ Qu terms. For the parameters associated with these 
terms, the difference between the results with and without the renormalization is not large. 
For instance, Xn of the Cr nuclei vanishes before the renormalization, and it remains small 
after the H^CM. In the = 32 nuclei, we have very small Xu both in the unrenormahzed 
and H^CM cases. A certain nucleus-dependence has been expected for the x parameters [Q. 
This variation with the increasing valence nucleon number is rapid in this region, because 
the size of the shell is small, compared with heavier nuclei. 

It is found that, while the Majorana interaction is negligibly small before the renor- 
malization, it becomes sizably repulsive due to the renormalization. In reality, the totally 
symmetric states are pushed down more than the mixed-symmetry states, absorbing the 
more effect of non-SD degrees of freedom. Thereby the mixed-symmetry states are pushed 
up to some extent, relative to the lowest-lying states. This point is already viewed in the 
fermion spectra shown in the preceding section. The repulsive Mi term after the H^CM is 
compatible with increase of the excitation energy of the l"*" states in Figs.|I], 0, ^ and 



It has been shown ||8|j20|l that, if we calculate a boson image of a schematic interaction, 
the Majorana terms emerge as a renormalization effect. A similar effect occurs also for 
a realistic interaction derived from the G-matrix. Nevertheless, there is a difference in 
mechanism. Whereas only the influence of g'-boson degree of freedom is considered in Refs. 
PPUj], various correlation is included in the present renormalization. In particular, the core 



excitation effect looks to play a certain role in the present case. Energies of the lowest O"*" 
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and 2"*" {i.e., symmetric) states are lowered greatly, mainly due to the coupling to the core 
excitation. This mechanism works less for the mixed-symmetry components, resulting in 
the repulsive Majorana terms. 

Comparing the H^CM results on the Majorana terms among the four nuclei, we find 
certain nucleus-dependence of the parameters. In ^^Fe, is fairly large with ^2 and ^3 
remaining quite small. The other nuclei have positive values for all of these parameters. 
The ^2 parameter is somewhat smaller than the others. 

It is also noticed that and the last term of Eq.(|2^) are not negligibly small. Though 
they hardly contribute to the lowest-lying states in these nuclei, some higher-lying states 
are affected to a certain extent. 

By diagonalizing the IBM-2 hamiltonian, we obtain the energy levels within the IBM-2. 
They are already displayed in Figs.|l|, 0, |] and |, by using the parameters after the H^CM. 
Since we include all the one- and two-body terms in the boson hamiltonian, the IBM-2 
levels in ^^Fe, where A^^ = 2, are exactly the same as those of the collective fermion space. 
In ^^Cr, the boson energy levels are very close to those of the collective fermion space, while 
the perfect agreement can be made if we use three-body terms in the boson hamiltonian. 
This indicates that the boson many-body terms are not important. The same holds for ^^Fe 
and ^^Cr. 

We next turn to electromagnetic transition operators. The following shell-model E2 
operator is assumed, 

T{E2)= Y: efj:r^Y^'\r,)^ (28) 

with ef = 1.4e and ef = 0.9e. The single-particle matrix-elements are evaluated by the 
harmonic-oscillator wavefunctions with b = 56^/^ = 1.956fm. For the Ml operator. 



T(M1) = 




and we take gf^ = Qi p^, g^p = O.^gg p'^. These electromagnetic operators are the same as in 
Ref. m. While medium effect on the M3 operator has not been explored sufficiently, there 
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is an evidence in the s(i-shell that no quenching is necessary to describe M 3 transitions ^1 . 
The shell-model M3 operator is taken to be equal to the bare-nucleon operator, 



p=TT,iy [ iep iep ) 

where the (^-parameter are the same as for the Ml operator. 
The IBM-2 operators are 

T^{E2) = 5: {[dls, + .tJjW + x'Md.f^} , (31) 

p=TT,U 

T'^m) = \/J E Jf = Vw[dld,f\ (32) 

V p=-n,v 

T^iMS) = E Mdld.f^- (33) 

p=1T,U 

The parameters introduced above are evaluated from the matrix elements within the collec- 
tive fermion space, analogously to the mapping for the hamiltonian. The resultant lBM-2 
parameters are shown in Table |ITT| . 

The renormalization enhances the boson effective-charges (ejf), gaining more quadrupole 
collectivity. In phenomenological studies, = is sometimes postulated. This is sup- 
ported by the present microscopic study. Though the shell-model effective-charge is smaller 
for neutrons than for protons, the neutrons have more quadrupole collectivity because the 
size of the valence shell is bigger. As a consequence, and are not so different. There 
is no apparent reason for the x' parameters to agree with the x parameters which appear in 
the boson hamiltonian, because we adopt a realistic shell-model interaction, not a schematic 
proton-neutron interaction like ■ Qu- Nevertheless, x' is close to x iii any case. It has been 
expected that ~ gj'^^^^ = 1 and g^ ~ gj'^J;'^ = 0, since the nucleon-spin degrees of freedom 
are not so active in the quadrupole collective states. However, g^ is certainly larger than 
unity in the result without the renormalization. This happens because the proton S'D-pairs 
consist of the single-j orbit of 0/7/2- On the contrary, when we carry out the renormalization, 
the excitation from O/7/2 to O/5/2 leads towards saturation of the nucleon spin. Therefore 
the g^ value is reduced until ~ 1 is restored to a good extent. The quenching of P^^j, due 
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to the renormalization is explained by a similar spin-saturation mechanism. It is somewhat 
surprising that the influence of the renormalization on and (3s^i, is so small. This will be 
partly because the spin saturation occurs already in the unrenormalized neutron parameters, 
resulting in almost vanishing values. 

The same quenching mechanism for the magnetic transition parameters will prevail in 
heavier nuclei, where a unique-parity orbit is involved in the valence shell. Since the spin- 
orbit partner of the unique-parity orbit is absent in the valence shell, the nucleon-spin content 
could influence and (3s^p to a considerable extent, when we ignore effects of excitation 
across major shells. However, by taking into account the excitation effects, those parameters 
will be quenched, owing to the spin-saturation tendency in the dynamics. It is expected 
for the Ml parameters that approaches unity and g^ almost vanishes. 

The dependence of the IBM-2 parameters on valence nucleon numbers is accounted for, 
in most cases, in terms of the quasi-spin properties of the relevant nucleon operator [Q. The 
E2 operator behaves as a vector in the quasi-spin space . Weak dependence is suggested 
for e^, and is confirmed in Table |ITT[ The number-dependence of x' is strong, as is expected. 
So is X ill the boson hamiltonian. Since T(M1) and T(M3) are quasi-spin scalars [^], the 
g^ and /S^^p parameters are expected to be nearly constant. This is true for g^, g^ and /Ss^-,,, 
but a considerable deviation is seen in Ps^^. This is a sort of many-body effect, originating 
in the subshell structure. 



VI. DISCUSSION ON H"CM 
A. Choice of primary bases 

In this section, we return to discussion on the H"CM. 

Even when the subspace wf'^ is fixed, there still remains an ambiguity in choosing the 
primary orthonormal bases }. The bases can be changed by a unitary 

transformation. Since the couplings between (p^^ and (f)^^^ are not uniform for various 
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A, the unitary transformation may lead to a different renormalized basis-set in the H"CM. 
Moreover, because of the orthogonahzation stated in Section \[T\\ , the H"CM bases generally 
depends on the ordering of the primary bases. 
The following three choices will be possible: 

(i) Some orthonormal basis-set is postulated for by a physical insight. 

(ii) The eigenstates within Wj'^ are taken as "^f^- They are placed in order according to 
the eigenenergies. 

(iii) \E'^°'' is redefined for each step n, so that ¥^ should be an eigenstate within wj" ^\ 
They are put in order according to the eigenenergies. 

To fulfill the space closure of Eq.(^, the couplings between renormalized states and re- 
maining degrees of freedom should be small. From this viewpoint, the latter choice seems 
favorable. On the other hand, the latter requires more complication in the numerical treat- 
ment. 

The basic dynamical properties should be well represented by the primary bases, other- 
wise the renormalization does not converge with small n. In the practical case of the Cr-Fe 



nuclei in Sections |T| and |V|, we have adopted (|), with the UTr+ui^) ® SUf{2) bases of 



IBM-2 |2J]. This will be appropriate because those nuclei seem to be nearly spherical. In 



deformed region, another choice might be better 



B. Choice of renormalized bases 

One could pursue the convergence of H"CM, by increasing the power n. However, besides 
tedious computations, it might lead to a dissimilar wavefunction of from that of '^f^- 
Then it is not reasonable to regard as a renormalized state. In some cases, this is 
circumvented if an eigenstate in F^"^ having the largest overlap with is adopted as 
instead of the lowest-lying one. 
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A complication occurs when there is a substantial fragmentation of the state which 
carries the main character of the primary state. As has been pointed out in Ref. [|I|, the 
mixed-symmetry 2+ state of ^^Fe may be the case. The 2^ and 2^ states share appreciable 
fractions of this collective component. As n increases, the diagonalization within F^"^ will 
lead to a problem at some value of n; plural eigenstates have considerable amplitudes of the 
primary basis ^1'^'^'*. It will not be desirable, in such a case, to choose a single eigenstate in 
F^"^ as a renormalized basis ^I/^"^. 

A solution to this problem is to adopt a linear combination of the few eigenstates in F^""^ 
as a renormalized basis. We can set a criterion of minimum amplitude for the states to be 
included. Another practical choice is just stopping at a certain n. It should be emphasized 
that, in any case, monitoring the H"CM outcome for each step is significant. 

In the actual case of the Cr-Fe nuclei, we do not come across the problems stated above, 
up to the H^CM. The lowest-lying eigenstate in F^"^ has the largest overlap with \ and 
no serious fragmentation is viewed. As shown in Section 0, the convergence in 0^ and 2^ 
is so rapid that we could acquire a good approximation by the H^CM. In this respect, the 
H^CM seems good enough to investigate collective states of the Cr-Fe nuclei. 



C. H"CM and Lanczos method 

As has been mentioned earlier, there is a common part between the H"CM and the 
Lanczos diagonalization method. 

The H"^CM energy levels are obtained via two steps of diagonalization; one within the 
subspace F^"^ and the other within Wj^\ The dimension of F^"^ is (n + 1), while that of 
Wj^'^ is /. The basis production and the diagonalization within f1^"^ is similar to the Lanczos 
method. In the case that there is a single basis in wf \ the H'^CM procedure is the same 
as the Lanczos method starting from since we do not need the diagonalization within 
Wf\ In other cases, the orthogonalization between F^"^ and F^^ (A 7^ A') in the H"CM 
does not appear in the Lanczos method. Apart from this difference, the H"CM is exploited 
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so as to make good use of the advantage of the Lanczos method. 

An emphasis should be put on the primary bases: we have requested that they should 
have a simple structure but carry the basic dynamics of the system. For instance, the SD 
states are taken as the primary bases in the application to the Cr-Fe nuclei in Section 0. 
Because of these properties, it is expected in many cases that the H"CM is more efficient 
than the Lanczos method. The number of bases in the H"CM is given by l{n + 1) for each 
J, which can be smaller than that necessary in the Lanczos method. In practice, within a 
fixed J, even less than 10 bases yield good accuracy for the lowest-lying levels of the Cr-Fe 
nuclei via the H"CM, whereas about 50 bases are normally required in the Lanczos method. 
It is noted that each diagonalization is performed for a matrix with quite a small dimension, 
(n + 1) or /. 



D. H"CM and a perturbative renormalization 

In this subsection, the H'^CM is discussed in connection with a perturbative method of 
renormalization. A more detailed discussion is given in Ref. f^. 

Let us recall Feshbach's projection method [|lOl, a well-known method of incorporating 



truncation effects. We define the P-space as the space of the primary bases; the original 
SD space in the present case. The total space corresponds to the k < 2 shell-model space, 
while the original hamiltonian is the Kuo-Brown hamiltonian. We shall take into account 
the effects of the Q-space, which is spanned by the non-primary bases. Note that all the <p- 
bases of Eq.(|) belong to the Q-space. The projection operator onto the P-space is denoted 
by P, and that onto the Q-space by Q, namely Q = 1 — P. 

In Feshbach's method, the renormalized P-space hamiltonian is given by 

Hp = Hp + PHQ—^-—QHP, (34) 

E - Hq 



where 



Hp = PHP, Hq = QHQ. (35) 
26 



For the exact treatment of Hq, eigenvalues of Hq have to be calculated. Let an eigenstate 
of Hq be denoted by \qi). Notice that Q\qi) = \qi) and Q = J2i lli) {lil- Then the second 
term on the right-hand side of Eq.(0) is rewritten as 

1 



(36) 



E - E{qi 

where E{qi) = {qi\H\qi). Though Eq.(0) gives an exact way to incorporate the influence 
of the Q-space into the hamiltonian, it is difficult and not advantageous to handle without 
any approximation in most cases, because of the following two reasons. If one wish to know 
exact eigenenergies, Hq must be treated exactly, which is usually a matrix with enormous 
dimension. Moreover, a non-linear coupled equation must be solved, since the eigenenergy 
E is also contained in the denominator of the second term. 

For the sake of simplicity, our discussion is restricted to a fixed J (conserved quantum 



number), without loss of generality. We introduce the following state generated from \E' 



(0) 

A 5 



\) = H^^^ki) oc QH 



, (A = l,2,---,0 



(37) 



where 



(0) 



HQH 



(0) 



(38) 



It is noticed that, with the notation in Section |T|, (j)^ can be expressed as (o) (o) (o), • 
H'$^^\ which is different from the basis (p^^'^ only in the lack of orthogonalization between 
and '^y (A ^ A'). 

We now substitute a c-number E{(f)^ = {(f))^\H\(f)^) for Hq in Eq.(p^. Then a diagonal 
matrix element of Hp is approximated by 

1 



(0) 



H, 



(0) 



E - E{4)^) 

,(0)| WhT,(0)\ 



{^''■;;'\H\^>';;>) + 



,(0)|p,- 

1 



^i¥^''\HQH\^>T)- 



E- E' 



(39) 



This is a kind of closure approximation, since it is given by replacing the energy denominator 
by a c-number. It is not easy, in general, to evaluate E in the energy denominator properly. 
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By substituting unperturbed energy E 



(0) 



(0) 



H 



"^^x^) for it, Eg ■ (|39D becomes equivalent 



to the second-order perturbation combined with the cfosure approximation. It is remarked 
that, however, the unperturbed energy is too high in the practical case of the Cr-Fe nuclei, 
causing too big renormalization effect. If we neglect the non-orthogonality between (j)^ 
and 0y for A 7^ A', the space fI^^^ becomes |\E'^°'', In addition, if E is estimated by 
diagonalizing the hamiltonian in this F^^'' space, we obtain E^^^ . Then the right-hand side 
of Eq.(|5^) is equivalent to the diagonal element of the H^CM collective hamiltonian. Note 



that E^^ is lower than Ef'\ and is closer to the exact E for low- lying states. 

Observing the above relation between the H^CM and the perturbative renormalization 
with closure approximation, we can claim that the H^CM is an improvement from the 
perturbative method in the following points; (i) the over-counting arising from the non- 
orthogonality between (p^ and (py {X ^ A') is removed, (ii) off-diagonal elements are evaluated 
in a consistent manner with diagonal ones, and (iii) E in the energy denominator is improved. 
If the overlap between cp^ and (py is negligible, which somewhat depends on how to choose 

and (p^x^ becomes quite similar. In order to 



the original basis-set j^E'^'''; A = 1,2, 



further make the relationship between the two methods more transparent, we shall restrict 

ourselves to the diagonal elements and ignore how E is estimated. 

If the second term of the right-hand side of Eq. ( p4D is expanded by the parameter 

EM - E{q,) 



i.X 



E - Ei<p, 



(40) 



we obtain 



(0) 



PHQ 



E-Hc 



-QHP 



E-E( 



AO) 



HQH 



H 2 



Q) 



E-Ei 



where 



\Hl\cPx) - [E{cPx)f 



(41) 



(42) 



represents variance of Hq in the state (py In a similar manner, the 0(C") term corresponds 
to correction due to the n-th moment of the distribution of |q'j)'s. By comparing Eq.(PD 
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with Eq.(^), it is found in O(C^) that the distribution of (l)^ o'^^r the eigenstates of Hq 
generally enhances the effect of the renormalization on energies. 

The above discussion is useful to acquire an intuitive picture of the H"CM. The terms 
regarding Hq are fully taken into account in the H^CM. Therefore, as far as diagonal elements 
are concerned, the difference between Eq.(0) and the H^CM is only in 0(C^)- In comparison 
with the H^CM, an advantage of the H^CM is the inclusion of the O(C^) effect [Q. As stated 
already, in the perturbative theory, which is connected to the H^CM well, the coupling of a 
primary state with the outer space (Q-space) is treated by using an averaged energy of the 
non-primary states E{(f)x), ignoring their distribution. The second-order effect is, in essence, 
the correction due to the distribution of the coupled states in terms of the variance, as shown 
in Eg . (Pl) . By extending the present discussion to higher order, it turns out that the n-th 
order effect of H"CM essentially corresponds to the n-th moment of the distribution of (p^^ 
in the Q-space. 

We now look back at the results shown in Figs.|l], |], |^ and ^ in Section In proceeding 
from H^CM to H^CM, the higher-lying states tend to go down more sharply. This implies 
that the variance axiHq) is more important in the higher- lying state. On the other hand, 
some mixed-symmetry states with relatively low energy (for instance, the lowest 1"*" and 
the second 2+ states in the collective space of each nucleus) do not come down so rapidly, 
compared with their surrounding states. Relatively small ax^Hq) is suggested for those 
mixed-symmetry degrees of freedom, giving rise to the repulsive Majorana interaction in the 
IBM-2 hamiltonian. 

VII. SUMMARY 

In order to study quadrupole collective modes based on a realistic shell model, we de- 
velop the if"-cooling method (H"CM), which leads to a wavefunction renormalization by 
incorporating the effect of the dynamical correlation. In practice, the H"CM is applied to 
the SD space of the Cr-Fe nuclei; ^^Fe, ^^Cr, ^^Fe and ^^Cr. While the shell-model ground- 
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state wavefunction is not fully covered with the simple 5'D-pair degrees of freedom, the 
shell-model O^" and 2f energies and wavefunctions are nicely approximated by considering 
up to the second power of H (H^CM). On the other hand, as far as the energy difference 
is concerned, the excitation spectra after the H^CM do not differ very much from those 
without the renormalization. Note that the ^^Ni-core excitation is taken into account in the 
present renormalization, as well as some effect of other like-nucleon pairs. 

An extended OAI mapping is also developed and apphed to the Cr-Fe nuclei. This is the 
first work to evaluate the IBM-2 parameters from a realistic shell-model hamiltonian. The 
wavefunction renormalization is converted to a renormalization of the IBM-2 parameters. 
Some effects of the renormahzation are discussed. Although most parameters in the IBM-2 
hamiltonian do not change considerably, the Majorana interaction becomes sizably repulsive 
as a renormalization effect. It is indicated that many-body terms are unnecessary in the 
IBM-2 hamiltonian. In the transition operators, the H"CM gives rise to spin quenching for 
the Ml and M3 proton parameters, as well as E2 effective-charge enhancement. The x 
parameters in the E2 operator are shown to take close values to those in the hamiltonian. 
This situation is not influenced by the renormalization. 
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TABLES 

TABLE I. Overlaps of wavefunctions of lowest-lying collective states, before and after the 
renormalization via the H'^CM, with those of the shell-model eigenstates (%). 



Nucleus 


unren. 


0^ 

J. 


H^CM 


unren. 


25^ 

1 


H^CM 


56 pg 


55.6 




97.0 


49.8 




93.3 


54Cr 


53.5 




95.8 


47.5 




92.3 


58Fe 


46.0 




91.2 


36.8 




83.6 


56Cr 


54.6 




95.1 


47.4 




90.9 
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TABLE II. Parameters for IBM-2 hamiltonian derived from the Kuo-Brown shell-model hamil- 



toiiiau. 



Parameter 


56 


Fe 


54 


Cr 


58 


Fe 


56 


Cr 






unren. 


H^CM 


unren. 


H^CM 


unren. 


H^CM 


unren. 


H^CM 




(MeV) 


1.022 


1.362 


1.022 


1.671 


1.022 


1.216 


1.022 


1.591 




(MeV) 


1.426 


1.441 


1.482 


1.270 


1.888 


2.004 


2.178 


2.143 


K 


(MeV) 


0.839 


0.926 


0.679 


0.765 


0.696 


0.710 


0.558 


0.538 


Xtt 




—0.933 


—1.220 


0.000 


—0.202 


—0.933 


—1.194 


0.000 


—0.164 


Xu 




— 1.250 


— 1.099 


— 1.239 


— 1.150 


—0.005 


0.126 


0.013 


0.205 




(MeV) 


—0.065 


0.303 


0.084 


0.186 


—0.043 


0.135 


0.080 


0.192 




(MeV) 


0.000 


-0.021 


0.000 


0.127 


0.000 


0.095 


0.000 


0.098 




(MeV) 


-0.009 


0.061 


0.008 


0.171 


-0.031 


0.198 


-0.020 


0.152 


bo 


(MeV) 


0.330 


0.575 


0.123 


—0.099 


0.565 


0.212 


0.337 


-0.119 




(MeV) 


0.101 


0.320 


0.034 


0.100 


0.006 


-0.076 


-0.001 


-0.187 


h 


(McV) 


-0.064 


-0.306 


-0.104 


-0.130 


0.003 


-0.197 


-0.042 


-0.030 




(MeV) 






0.000 


0.012 






0.000 


0.060 


C2,7r 


(MeV) 






0.105 


-0.155 






0.100 


-0.020 


C4,7r 


(McV j 






-0.424 


-1.215 






-0.469 


-1.033 




(MeV) 










0.082 


0.063 


0.141 


0.233 


V2,v 


(MeV) 










-0.441 


-0.287 


-0.615 


-0.690 


CQ,v 


(MeV) 










2.066 


1.609 


3.399 


2.439 


C2,v 


(MeV) 










-0.308 


-1.347 


0.053 


-0.924 




(MeV) 










0.447 


0.170 


0.694 


0.593 
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TABLE III. Parameters for IBM-2 transition operators derived from the shell-model operators. 



Parameter 


56pg 

unren. H^CM 


s^Cr 
unren. H^CM 


58Fg 

unren. H^CM 


56 Cr 
unren. H^CM 


(efm^) 
(efm2) 


6.635 7.958 
6.034 7.273 
-0.933 -1.127 
-1.226 -1.307 


5.418 6.563 
6.009 6.618 
0.000 -0.336 
-1.205 -1.334 


6.635 7.609 
5.057 6.037 
-0.993 -1.178 
0.011 -0.049 


5.418 6.243 
5.004 5.365 
0.000 -0.331 
0.035 0.010 


y-K \H'N) 


-0.027 -0.037 


-0.039 -0.080 


1 1111 
-0.035 -0.015 


1 9iifi 1 1 in 
-0.038 -0.042 


P3,u (/^ivfm^) 


69.1 51.5 

-18.8 -20.7 


69.1 58.4 

-22.5 -26.3 


69.1 51.4 
6.1 8.6 


69.1 59.0 
5.3 10.1 
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FIGURES 

FIG. 1. Energy levels in the collective space (without renormalization, with renormalization 
via the H^CM, and with renormalization via the H'^CM), in comparison with the shell- model 
ones m ^^Fe: (a) Lowest 0+ and 2+ energy eigenvalues in each space, relative to the shell-model 
ground-state energy, (b) Energies relative to the lowest 0"*" level in each space. 

FIG. 2. Energy levels in the collective space compared with the shell-model ones in ^^Cr. 
Energy levels obtained from the IBM-2 hamiltonian are also shown. 



FIG. 3. Energy levels of ^^Fe. The experimental data are taken from Ref. [16|. The calculated 



energy levels are obtained by the k < 2 shell-model calculation with the Kuo-Brown hamiltonian. 



FIG. 4. Energy levels of ^^Cr. The experimental data are taken from Ref. |17|. The calculated 



energy levels are obtained by the k < 2 shell-model calculation with the Kuo-Brown hamiltonian. 

FIG. 5. Energy levels in the collective space, as well as the IBM-2 energy levels, compared 
with the shell-model ones in ^®Fe. 

FIG. 6. Energy levels in the collective space, as well as the IBM-2 energy levels, compared 
with the shell- model ones in ^^Cr. 
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